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Abstract 

We note that it is possible to construct a bond vertex model that displays g-state Potts criticality 
on an ensemble of 4> 3 random graphs of arbitrary topology, which we denote as "thin" random graphs 
in contrast to the fat graphs of the planar diagram expansion. 

Since the four vertex model in question also serves to describe the critical behaviour of the Ising 
model in field, the formulation reveals an isomorphism between the Potts and Ising models on thin 
random graphs. On planar graphs a similar correspondence is present only for q = 1, the value 
associated with percolation. 



1 Introduction 



A description of the Potts model on planar random graphs as an ice type vertex model on the associated 
medial graphs was developed some years ago by Baxter et.al. Q. In this paper we formulate a bond 
vertex description of the Potts model on non-planar random graphs, which we call "thin" graphs to 
distinguish them from the fat graphs which appear in a planar graph expansion. The study of spin and 
vertex models on such thin random graphs is of interest as a way of obtaining mean field theory results 
@! H' H' H without the boundary problems of the Bethe lattice or the inconvenience of infinite range 
interactions. Mean field behaviour is present because the thin graphs look locally like a Bethe lattice, but 
all the branches of the tree-like Bethe lattice structure are closed by predominantly large loops. Planar 
graphs on the other hand have a loop distribution that is strongly peaked on short loops and a fractal-like 
baby universe structure ||. 

Nonetheless, the methods used for calculating the partition functions of spin models on thin graphs 
may still be borrowed from the study of the planar random graphs . These methods are based on the 
observation that planar graphs can be thought of as arising as Feynman diagrams in the perturbative 
expansion of integrals over N x N Hermitian matrices. Each edge in such a Feynman diagram is "fat" 
or ribbon-like, since it carries two matrix indices. In the N — > oo limit fat graphs of planar topology 
are picked out, since there is an overall factor of N x for any graph, where x is the Euler characteristic. 
The N — > 1 limit, on the other hand, weights all topologies equally. In this case the matrix fat graph 
propagators degenerate to scalars, so we denote the generic random graphs of the N — ► 1 limit as "thin" 
graphs. The Feynman diagram approach used in studying the statistical mechanics of models on planar 
random graphs [] still applies for the thin graph case too. In fact, life is even easier in this case since one 
is dealing with scalar rather than matrix integrals. 

The partition function for a statistical mechanical model on an ensemble of thin random graphs with 
2m vertices can be defined by means of an integral of the general form g 

where the contour integral over the vertex coupling A picks out graphs with 2m vertices, S is an appro- 
priate action, K is the inverse of the quadratic form in this action, and the <pi are sufficient variables 
to describe the matter in the theory. The factor N m gives the number of undecorated (without matter) 
graphs in the class of interest and disentangles this usually factorial growth from any phase transitions. 
For the class of thin <fi 3 (3-regular) random graphs we discuss here 

/l^ (6m -1)11 

In the large m, thermodynamic, limit saddle point methods may be used to evaluate equ.(|j]). The saddle 
point equation for A may be decoupled by scaling it out of the action as an overall factor, leaving any 
critical behaviour residing in the saddle point equations for the matter fields (j>i . Phase transitions appear 
as an exchange of dominant saddle point. 



2 The Potts Model on Thin Graphs 

The Hamiltonian for a g-state Potts model can be written 

<ij> 

where the spins o"j take on q values. The action which generates the correct Boltzmann weights for the 
g-state Potts model on </> 3 graphs to be used in equ.(|l|) is || |8| 

s=^E^-cE^-^I> 3 ( 4 ) 

i— 1 i<.j 2—1 



1 Which is equivalent to coupling the models to 2D quantum gravity. 
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where c is 1/ (exp(2/3) + q — 2). Since for any q one finds a high temperature, disordered phase, solution of 
the form fa = 1— (q— l)c, Vi bifurcating to a broken symmetry, ordered phase, solution fa = . . . i 7^ 0i 
at c = 1/(2(7 — 1) an effective action with only two variables 4>,4>, where <p — fa, <f> — fa, fa, . . . ,fa, 
suffices to describe the transition || 

Si = \(q - 1) [1 - c(q - 2)] fa 2 ~ j{q ~ 1)0 3 + \tf ~ ~0 3 ~ c(q - 1)#. (5) 

In the high temperature phase </> = <j> and this collapses to 

S 2 = q -(l-c(q-l))fa-^. (6) 

The bifurcation point is not the first order transition point displayed by the model for q > 2, but 
rather a spinodal point. The first order transition is pinpointed by observing that the free energy in the 
saddle point approximation is the logarithm of the action S so the first order transition point is given by 
the c value, and hence temperature satisfying Si = S 2 - This gives the critical value of c as 

cid^pL, (7) 

q z 

in agreement with other mean field approaches || . 

The nature of the phase diagram can best be clarified by examining a diagram of the magnetisation 
defined by 

m = jz r (8) 

(f + (q- 1)03 j 

vs c, as plotted in Fig.l for q = 4. We can see that the bifurcation point at P, c = l/(2q — 1) for general 
q, lies below the first order transition point at Q, c = [1 — (q — l)~ 1//3 ]/(q — 2). The second spinodal 
point at O, c = [q—1 — 2^/q — T\/[(q —l)(q — 5)], is given by the vanishing of a square root in the saddle 
point solution. As q — * 2 O, P, Q coalesce, the skewed pitchfork of Fig.l becomes symmetrical and we 
recover the continuous mean field transition of the Ising model - which is equivalent to the q = 2 Potts 
model. This is to be expected since setting q = 2 directly in the action of equ.(||) and the magnetisation 
of equ.(^) recovers the Ising action and magnetisation. 

Since q appears explicitly as a parameter in equ. (J5J) the formalism is ideally suited to studying perco- 
lation on thin graphs as the q — * 1 limit of the Potts model, as well as the random resistor, q — > 0, and 
dilute spin glass, q — >• 1/2, problems. 



3 A Potts Vertex Model 



In order to recast the Potts action as a bond vertex model we first carry out the following rescaling on 
equ. (pi) 
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S = 1(X 2 + Y 2 ) - [X 3 + 3 K *XY 2 ] 



y/(q-l)(l-c(q-2)) 

cut (jf 
ing of 

A(k*) 3/2 



(9) 



which gives both the quadratic terms the canonical coefficient of 1/2, followed by the linear transformation 
<fi — > (X — Y)/\/2, 4> — > (X + F)/v2 and a further rescaling of the quadratic terms and A to obtain 



3 2 

where «;* = (! — «)/(! + k) and 



Y A + —YX Z 

K* 



(10) 



(g-l)V2(l - (q-2)cf/ 2 



1 c 2 (q-l) 
l-(q-2)c 



(11) 
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The notation in equ.(10) has been chosen to facilitate comparison with the Ising model in held below. 

We can see from the above that the action of equ.(|^) for the Potts model on thin 3 graphs can be 
transformed via some rescalings and a linear transformation of the variables into a 4- vertex model on the 
4> 3 graphs. What is more, the vertex model is symmetric since the weight depends only on the number 
of X and Y bonds at a vertex The different vertices in the model are shown in Fig. 2. Although the 
result of a straightforward transformation, equ.(jl0|) entails a surprising consequence. The action for the 
Ising model (i.e. the q = 2 state Potts model) in held on thin graphs ^ is given by |n| 



Tr 



Y 2 ) - gXY + ^ [e h X 3 + C - ! 'Y 3 ] , 

o 



where g = exp(— 2(5) and h is the external field. Carrying out the transformations 



(12) 



X 
Y 



(X + Y)/V2 
(X-Y)/V2 



followed by the rescalings X 
four-vertex model 



X/(l-g)V*, Y -> 57(l + 0)Va A 



(13) 

\/2A(l — g) 3 / 2 again gives a 



Y 2 )- 



A cosh(/i) 



[X 3 + 3g*XY 2 ] - 



Asinh(/ l )(< 7 *) 3 / 2 



Y A + — X Z Y 

9* 



(14) 



where g* = (1 - g)/(l + g) 

We thus find that the vertex model action for the q state Potts model and that for the Ising model in 
field are isomorphic under the identifications 



tanh(/i) 
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1 + v 



(15) 



In fact, backtracking to equ.(||) we can see that the rescaling of equ.(||) transforms the action into that 
of an Ising model in field, even before the transformation to a vertex model 



5' 



+ 4?) 



Xv 



A 7q 



(16) 



The equivalence between the Potts and Ising models is surprising from the physical point of view 
since we know already that the Potts models display a first order transition for q > 2 || |j| whereas the 
Ising model displays a (mean-field) second order transition. Things become clearer when we consider the 
mapping of the Potts values for c(q, (3) at the critical and spinodal points onto the parameters k, v, which 
play the role of temperature and field in the Ising model (or more precisely exp(— 2/3) and exp(2/i)). The 
corresponding values of k for 0,P,Q are shown in Fig. 3 and for v in Fig. 4. We find that the hrst order 
transition point in the Potts models at Q, where c = [1 — (q — l) -1 / 3 ]/ (q — 2), maps onto 



= ± 
= 1 



(g _ 1)2/3 _ ( g _ 1)1/8 



(17) 



where the sign on the right hand side of the expression for k is chosen to give a positive answer depending 
on whether q <> 2. Remarkably we see that this corresponds to a zero-field point in the Ising model, 
with k <C k c = 1/3, the Ising critical value. Since this means that (3 > critical ising-, this point lies on 
the zero field line separating the two possible spin orientations in the ordered phase, so the first order 

2 On random graphs the notion of orientation is lost, so there are fewer possibilities for defining independent vertex 
weights than on regular lattices. On planar random graphs one can still define a cyclic ordering of bonds round a vertex 
consistently, but even this is lost on thin graphs. 

3 And on planar graphs too, if we take X, Y to be matrices. 
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temperature driven transition of the Potts model is mapped onto the field driven transition of the Ising 
model. As q — > 2 from above or below, k — > 1/3, the Ising critical value, and the transition becomes the 
continuous mean-field Ising transition. 

The spinodal point at P where c = l/(2q — 1) maps onto 



q-1 



(2g-l)(g+l) 
2q- 1 



(18) 



(q- l)V2(g+l) 

which tends towards the standard mean field Ising transition point in zero-field as q — > 2. Similarly the 
other spinodal point at at O, c = [q — 1 — 2-y/g — lJ/[(<7 — l)(q — 5)], maps onto 

(<z-i) 1/4 



(2^9 = 1 + 1)1/2(^/9=1 + 2)V2 
( v ^T+2) 3 / 2 (g-l) 1 / 4 
(2Vg=T+l) 3 / 2 



(19) 



so we see that both these points are not, in general, zero field points in the Ising transcription except at 
(7 = 2. 

A natural question to ask is what is so special about Q from the Ising point of view. If we look at the 
effect of the scaling of equ.(g) on the original Potts action of equ.(||) and demand a c value which gives 
v = 1 (i.e. zero field) we get precisely c = [1 — (q — l)~i/ 3 ]/(g — 2). The first order transition point at c( 
Q ) is thus the only point which maps onto zero field in the Ising model for general q. 



4 Discussion 

The main conclusion of this paper may be simply stated: the Ising model in field and the Potts model on 
thin </> 3 random graphs may be mapped onto one another. This can be seen either by mapping them both 
to a four vertex model, or by directly rescaling the Potts action. The equivalence is at root due to the 
fact that the Potts action of equ.(||) requires only two variables, just like the Ising model, even though 
the Potts spins have q states. We have also seen that the first order nature of the Potts transitions for 
q =/= 2 and the continuous transition of the Ising model are not in contradiction, since the Potts transition 
maps onto the field driven transition of the Ising model. The first order Potts transition point lies on the 
zero field Ising locus at j3 > ^critical ising and moves to the continuous Ising transition point as q — ► 2. 

One might enquire as to whether similar relations could exist on planar random graphs where the 
scalars in actions such as equ.(Q) are replaced by N x N Hermitian matrices. Unfortunately, a two variable 
(in this case matrix variable) effective action such as that in equ.(^J) does not appear to exist in such a 
case. However, one can arrive at a 3q+ 1 vertex model for planar graphs by using Kazakov's approach || 
of introducing an auxiliary matrix to decouple the Potts interactions. This transforms the matrix action 
for the q state Potts model 

( i—1 i<j i—1 J 

to 

{ i=l i 

which after a shift in fa and a rescaling may be written as 
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where v = exp(/3) — 1. In general this will not be equivalent to an Ising model since we have too many 
matrices, but when q = 1, which is related to the problem of percolation, we do recover an Ising- like 
action. 

It would be of some interest to see whether this transcription might shed some light on either perco- 
lation or the Ising model in field on planar random graphs. It would also be interesting to explore the 
relation between the planar graph Potts vertex model of equ. (|22|) and the medial graph vertex model of 
Baxter et.al. 
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Figure 1: The magnetisation m for a 4 state Potts model as calculated from the saddle point solutions. 
The high temperature solution is shown dotted, one low temperature solution solid and the other dashed. 
The first order jump in m is at Q. 
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Figure 3: k vs q for the points 0,P,Q. The curves touch at q = 2, k = 1/3 which is the Ising critical 
point. 
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Figure 4: v vs q for the points 0,P,Q. Since v = 1 corresponds to zero-field in the equivalent Ising 
model, we see that Q is a zero- field point for all q. For O, P on the other hand v = 1 only at q = 2. 
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